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Abstract 

Diagnosis and prognosis are necessary tasks for system re- 
configuration and fault-adaptive control in complex systems. 
Diagnosis consists of detection, isolation and identification of 
faults, while prognosis consists of prediction of the remain- 
ing useful life of systems. This paper presents a novel inte- 
grated framework for model-based distributed diagnosis and 
prognosis, where system decomposition is used to enable the 
diagnosis and prognosis tasks to be performed in a distributed 
way. We show how different submodels can be automati- 
cally constructed to solve the local diagnosis and prognosis 
problems. We illustrate our approach using a simulated four- 
wheeled rover for different fault scenarios. Our experiments 
show that our approach correctly performs distributed fault 
diagnosis and prognosis in an efficient and robust manner. 

1. Introduction 

Systems health monitoring is essential to guaranteeing the 
safe, efficient, and reliable operation of engineering sys- 
tems. Integrated systems health management methodologies 
include fault diagnosis and prognosis mechanisms, where di- 
agnosis involves detecting when a fault has occurred, isolat- 
ing the true fault, and identifying the true damage to the sys- 
tem; and prognosis involves predicting how much useful life 
remains in the different components, subsystems, or systems 
given the diagnosed fault conditions. The information on the 
fault size and its expected impact on system life can be used 
to initiate recovery and reconfiguration actions that mitigate 
the fault or extend system life. 
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A large body of research exists for both model-based di- 
agnosis (Gertler, 1998; Blanke et al., 2006) and prognosis 
methods (Luo et al., 2008; Saha & Goebel, 2009; Orchard 
& Vachtsevanos, 2009). However, the integration of diag- 
nosis and prognosis algorithms is seldom studied. In fact, 
many diagnosis methodologies leave out the fault identifica- 
tion step, which is necessary to perform a prediction from the 
current system state. Recently, we presented an integrated 
model-based framework for diagnosis and prognosis of com- 
plex systems, in which we made use of a common modeling 
framework for modeling both the nominal and faulty system 
behavior (Roychoudhury & Daigle, 2011). 

In (Roychoudhury & Daigle, 2011), the nominal system be- 
havior is estimated using an observer built with the nominal 
model. Faults are detected when a statistically significant de- 
viation between the nominal estimates and the observed mea- 
surements is observed (Biswas et al., 2003). Fault isolation 
compares the observed measurement deviations against pre- 
dictions of how the measurements would deviate for each 
possible fault (Mostennan & Biswas, 1999). Fault identifi- 
cation performs joint state-parameter estimation using multi- 
ple observers, where, for each fault, the faulty system model 
is constructed as the nominal model integrated with a hy- 
pothesized fault model (Roychoudhury, 2009). The prog- 
nosis module uses, for each fault hypothesis, a prediction 
model based on its faulty system model and the identified 
fault parameters, to predict the remaining useful life of the 
system (Daigle, Saha, & Goebel, 2012). However, this inte- 
grated solution performs the diagnosis and prognosis task in a 
centralized fashion, which is prone to single points of failure, 
and does not scale well as the size of the system increases. 

To overcome such problems, in this work, we leverage re- 
cent results for distributed diagnosis (Bregon et al., 2011) 
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and distributed prognosis (Daigle, Bregon, & Roychoud- 
hury, 2012), which make use of structural model decompo- 
sition techniques, to provide a systematic approach to dis- 
tributing the different diagnosis and prognosis steps presented 
in (Roychoudhury & Daigle, 2011). 

Distributed diagnosis is achieved by designing local dis- 
tributed subsystems based on global diagnosability analysis 
of the system, thus computing globally correct distributed 
diagnosis results without the use of a centralized coordina- 
tor (Bregon et al., 2011). These local distributed subsystems 
are then used to construct local event-based distributed diag- 
noses for distributed fault isolation. Distributed fault iden- 
tification is achieved by developing independent local state- 
parameter estimators for each hypothesized fault. Regarding 
distributed prediction, in (Daigle, Bregon, & Roychoudhury, 
2012) we developed an architecture that enables a large prog- 
nosis problem to be decomposed into several independent lo- 
cal subproblems from which local results can be merged into 
a global result. 

The main contribution of this paper is an integrated frame- 
work for distributed model-based diagnosis and prognosis of 
single faults based on structural model decomposition. The 
proposed framework scales well and the resulting subprob- 
lems are typically small and easy to solve, resulting in an ef- 
ficient and scalable distributed solution to the combined diag- 
nosis and prognosis problem. We perform a number of exper- 
iments on a simulated four-wheeled rover testbed (Balaban et 
al., 201 1) to demonstrate and evaluate our approach. 

The rest of the paper is organized as follows. Section 2 pro- 
vides the problem formulation for our diagnosis and progno- 
sis framework. Section 3 describes the distributed architec- 
ture and Section 4 briefly introduces its different components. 
Section 5 presents the case study and experimental results. 
Finally, Section 6 concludes the paper. 

2. Problem Formulation 

The nominal system model is represented as follows: 

x(t) = f(t,x(t),0(t),\l(t),v(t)), 
y (t) = h(t,x(t),9(t),u(t),n(t)), 

where x(t) G M nx is the state vector, 9(t) G R ne is the 
parameter vector, u (t) G M nu is the input vector, v(t) G 
is the process noise vector, f is the state equation, y (t) G M"* 
is the output vector, n(t) G R 71 " is the measurement noise 
vector, and h is the output equation. 1 

Faults in the system are represented as changes in the above 
nominal system model. In this work, we only consider sin- 
gle faults occurring as changes in system parameters, 0(t). 
We denote a fault, / G F, as a tuple, ( Q,gj ), where, 0 G 0 

'Here, we use bold typeface to denote vectors, and use n a to denote the 
length of a vector a. 


is the fault parameter, and g / denotes the fault progression 
function , which models the way fault / is manifested in pa- 
rameter 0 , i.e.. 

Hi) = ), U (<), m /(<)), 

where x f (t) = [x(t), 9{t)] T , 9 f (t) = [9(t)\{9(t)}, 

4>f(t)] T , 4> } (t) € K 71 */ is a vector of fault progression pa- 
rameters, and m/(£) G M nm / is a process noise vector asso- 
ciated with the fault progression function. 

To develop our integrated diagnosis and prognosis frame- 
work, the faulty system model for fault / = (0,gj) is con- 
structed from the nominal system model by including the pa- 
rameter as a state and augmenting the state equation by in- 
cluding the fault progression function, i.e., 

x/(f) = ff(t,Xf(t),9f(t),\i(t),v(t)), 
y(t) = h(t,x(t), 9(t), u(t), n(f)), 

where. 


f(t,x(t),0(f),u(t),v(£)) 


x(i) ' 

. 9f(t,Xf(t),9f(t),u(t),m(t)) _ 


so 

•QS 


The goal of diagnosis is to: (/) detect a change in some 6 G 0\ 
(ii) isolate, under the single fault assumption, the true fault 
/ G F, i.e., both the parameter 9 that has changed, and its 
fault progression function gj \ and (Hi) identify (i.e. estimate) 
the fault by computing p(x/(£), 6 f(t) |y(0 : t)), where y(0 : t) 
denotes all measurements observed up to time t. 

The goal of prognosis is to determine the end of (use- 
ful) life (EOL) of a system, and/or its remaining useful 
life (RUL). For a given fault, /, using the fault estimate, 
p(x.f(t),0 f(t)\y(0:t)), a probability distribution of EOL, 
p(EOLf(tp)\y(0:tp)), and/or RUL, p(RUL/(f P )|y(0:*p)) 
is computed at a given time point t P (Daigle, Saha, & 
Goebel, 2012). Since the prognosis problem is stochas- 
tic, EOL/RUL are random variables and we represent them 
by probability distributions. The acceptable behavior of 
the system is expressed through a set of n c constraints, 
Ceolj = {ci)T=i> w ^re c* : R nx f x R n$ f x E n,< — ► B 
maps a given point in the joint state-parameter space given 
the current inputs, (x/(f), 0/(t), u(£)), to the Boolean do- 
main B = [0,1], where Cj(x/(f), #/(£), u(i)) = 1 if the 
constraint is satisfied (Daigle, Saha, & Goebel, 2012). If 
Ci(xf(t), 0/(t), u(t)) = 0, then the constraint is not satis- 
fied, and the behavior of the system is deemed to be unac- 
ceptable. These individual constraints are combined into a 
single threshold function T E0 L f '■ ^ nx/ x x R nu — » B, 
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defined as 

2BOL/(x/(i),0/(t),u(t)) = 

1, 0 6 {cj(x/(t), u(t))}^! 

0, otherwise. 

So, EOL f may be defined as 
EOL f (t P ) = 

inf {t 6 R :t>tp and Teol, (x/(t), 0/(t), u(i)) = 1}, 

i.e., EOL is the earliest time point at which the threshold is 
reached. RUL is expressed given EOL as 

RUL/(* P ) — EOL/(^p) — tp. 

3. Distributed Architecture 

For a large system, both the diagnosis and prognosis prob- 
lems are correspondingly large. A centralized approach does 
not scale well, can be computationally expensive, and prone 
to single points of failure. Therefore, we propose to decom- 
pose the global integrated diagnosis and prognosis problem 
into independent local subproblems. In this work, we build 
on the ideas from structural model decomposition (Blanke et 
al., 2006; Pulido & Alonso-Gonzalez, 2004) to compute lo- 
cal independent subproblems, which may be solved in paral- 
lel, thus providing scalability and efficiency. Model decom- 
position is not a new concept, and several approaches have 
been developed for purposes of system identification, estima- 
tion, learning, and diagnosis (Staroswiecki & Declerck, 1989; 
Pulido & Alonso-Gonzalez, 2004; Williams & Millar, 1998). 
Structural model decomposition allows decomposing a global 
model into a set of local submodels for which local diagnosis 
and prognosis problems can be directly defined. The global 
model of the system, denoted as M, is defined as follows. 

Definition 1 (Model). The model of a system, M, is a tuple 
M = ( X , 0, U , y, C), where X is the set of state variables 
of x, 0 is the set of unknown parameters of 0 y U is the set of 
input variables of u, Y is the set of output variables of y, and 
C is the set of model constraints of f, h, and EOL constraints 
of Ceol /• 

The basic idea of the model decomposition problem is to de- 
compose the global system model into a set of submodels sat- 
isfying given constraints, such that each submodel contains 
sufficient analytical redundancy to generate fault hypotheses 
from observed measurement deviations. A submodel in our 
framework is defined as follows. 

Definition 2 (Submodel). A submodel Mi of a system model 
M = (X,0,f/,y,C) is a tuple M { = ©,,!/<, If, C<), 

where X* C X, ©< C 0, Ui C X U U U Y y and Y { C Y are 
the state, parameter, input, and output variables, respectively. 


and Ci C C are the submodel constraints. 2 

Next, we discuss the fundamental ideas of our model decom- 
position approach and show the constraints needed to obtain 
the different submodels for distributed diagnosis and progno- 
sis. Then, we propose our integrated approach. 

3.1. Model Decomposition for Distributed Diagnosis and 
Prognosis 

Model decomposition in our framework can be accomplished 
by using some variables (either measured variables or vari- 
ables for which the values are known) as local inputs, Ui , such 
that each one of the submodels satisfies a set of constraints 
and contains the minimum number of equations to compute 
a given set of outputs, Y{. Asa result, submodels computed 
this way contain only a small subset of the equations of the 
model that are decoupled from the rest of the system model 
equations. In general, any set of variables in the system can 
be chosen as the local inputs to the submodel, U{. The choice 
of Ui and the constraints to fulfill depends on the particular 
problem to be solved through model decomposition. 

The first model decomposition problem computes minimal 
submodels from the nominal system model. For this model 
decomposition problem, constraints are that submodels use 
the global model inputs and some measured values as local 
inputs, i.e., Ui C U U (Y - Yi). An algorithm for comput- 
ing the set of minimal submodels that satisfies these prop- 
erties is given in (Daigle et al., 2011 ), which is based on 
the model decomposition algorithms presented in (Pulido & 
Alonso-Gonzalez, 2004 ; Bregon et al., 2012 ). 

The second model decomposition problem computes sub- 
models for residual generation and fault isolation. For this 
model decomposition problem, the constraints are that sub- 
models are constructed by merging the minimal submodels, 
to fulfill global diagnosability conditions as in (Bregon et 
al., 201 1). 3 Once the globally diagnosable subsystems have 
been designed, the merged submodels are used for distributed 
residual generation and to compute event-based local diag- 
noses for fault isolation. These design and diagnoser com- 
putation processes are detailed in (Bregon et al., 2011 ). Be- 
cause the subsystems are designed to be globally diagnosable, 
the resulting local diagnosers are independent, and can pro- 
vide globally correct diagnosis results without a centralized 
coordinator. 

The third model decomposition problem computes submod- 
els for distributed fault identification. For each consistent 
fault hypothesis /, the joint state-parameter estimators are 
computed from the minimal submodels of the faulty system 
model M(f) with the constraints that Ui C UU(Y — Yi) with 

2 A faulty submodel for a fault / is denoted as A1i(/). 

3 In this work, a subsystem is globally diagnosable if all faults in the subsys- 
tem are distinguishable from every other fault in the system using only local 
measurements. 
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Distributed Diagnosis Distributed Prognosis 



Figure 1. An instantiation of the integrated diagnosis and prognosis architecture. 


Yi as a singleton. It will be shown later that the fault identi- 
fication module is the central part of our diagnosis-prognosis 
integration approach and provides the joint state-parameter 
estimations for the prediction module. 

Finally, for distributed prediction, the model decomposition 
problem starts off from the faulty system model, and, as de- 
tailed in (Daigle, Bregon, & Roychoudhury, 2012), it fulfills 
the following constraints: (i) the submodels use £/, C Up, 
where Up C XU U (here. Up is a set of variables whose 
future values can be predicted a priori , which depends on the 
hypothesized faults); and (if) each computed submodel has 
at least one c E CeoLj belonging to C*, and over all sub- 
models, all constraints in CeoLj are covered. This ensures 
that TeoLj may be computed for the system from the local 
constraints. 

3.2. Distributed Architecture 

Figure I illustrates an example architecture for our distributed 
diagnosis and prognosis scheme. At each discrete time step, 
k , the system takes as input both u k and y k and splits them 
into local inputs uj. and local outputs yj. for the local diag- 
noses. Within each M x local diagnoses nominal tracking is 
performed, computing estimates of nominal measurements, 
y k . The fault detector compares the estimated measurements 
against the observed measurements, to determine statistically 
significant deviations for the residual, = y\ — y k . Qual- 
itative values of the deviations in the residuals are used by 
the event-based diagnoser to isolate faults. The set of iso- 
lated fault candidates F^. together with the estimated nom- 
inal states, x^., parameters, 0 k> and the measurements, y*., 
are used as input for the corresponding identification mod- 
ule. Identification is performed for each hypothesized fault 
in a distributed way, e.g., for the isolated faults f\ and in 
Figure 1, we run an instantiation of the identification sub- 
model for each one the faults, i.e., M\(fi) and Mi(f2)- 
Fault identification uses the minimal submodels from the 
faulty system model, and computes local state-parameter es- 


timates p(xy 0 l j k \y l 0 . k )- These local estimates are then 
used as input to the prediction submodels. In some cases, 
the local estimates have to be split or merged with other 
estimates according to the prediction submodels. For ex- 
ample, in the figure, estimates from M\(f\ ), are used by 
both local prediction submodels and M^(fi), and 

those submodels may also need estimates not included within 
submodel M\(f\). These estimates are typically obtained 
from the local diagnosers or other fault identification blocks. 4 
Distributed prediction modules compute, for each hypothe- 
sized fault, local EOL/RUL predictions, p(EOL l j kp \y l Q :kp ) 
and v{RUL l j kp \y\ ) . kp ) y at given prediction time kp based 
on the local EOL constraints. Finally, local predictions 
are combined into global predictions p(EOLf, kp \yo :kp ) and 
p(RULf' kp \yo :kp ) for each hypothesized fault. The next 
section describes the details of the different modules of the 
distributed integrated diagnosis and prognosis architecture. 

4. Diagnosis and Prognosis Approach 

Figure 1 shows the basic modules of our distributed integrated 
approach. In this section we give details on how each mod- 
ule is implemented, and establish the integration between the 
diagnosis and prognosis tasks. 

4.1. Distributed Diagnosis 

For distributed diagnosis, each local diagnoser first takes a 
subset of the local inputs u k and local outputs yj., to compute 
an estimate of its output measurements yj.. Tracking is per- 
formed in discrete time using a robust filtering scheme, e.g., 
the extended or unscented Kalman filter (Julier & Uhlmann, 
2004), which provides accurate tracking in the presence of 
sensor noise, process noise, and discretization error. 


4 Si nee prediction submodels are constructed by using any variable which 
value can be hypothesized as input, in some cases, prediction submodels 
cannot always be formed by just merging the minimal estimation submod- 
els. To indicate this, we named our prediction submodels differently from 
the estimation submodels, e.g., M±(f \ ) and M^(}\) instead of M\(f\). 
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For fault detection, a statistical test is used to look for signif- 
icant deviations in the residual signal r k , which is computed 
as the difference between y k and the y l k . In our approach, we 
use a Z - test as described in (Biswas et al., 2003). 

Fault isolation is performed using local event-based di- 
agnoses, constructed as detailed in the previous sec- 
tion (Bregon et al., 201 1; Daigle et al., 2009). Fault isolation 
is triggered when a fault is detected, and it works as follows. 
Initially, all event-based local diagnosers start in their initial 
state, and the set of faulty candidates is empty. Local residual 
deviations cause the local diagnosers to move from one state 
to another. These residual deviations are abstracted to a tu- 
ple of qualitative symbols (<7i,cr 2 ) for each residual signal, 
where a\ represents magnitude changes and <r 2 represents 
slope changes. A -1- (resp. — ) value indicates a change above 
(resp. below) normal for a measurement residual or a positive 
(resp. negative) residual slope. A 0 implies no change in the 
measurement value or a flat residual slope. The symbols are 
generated using a sliding window technique as described in 
detail in (Biswas et al., 2003). If there is a match between 
an event from the current state and a tuple of qualitative sym- 
bols generated by any residual, the local diagnoser moves to 
the next state and remains active. If not, the local diagnoser 
blocks. This process continues until a local diagnoser reaches 
an accepting state, which corresponds to a unique isolation 
result. 

In our distributed diagnosis approach, identification submod- 
els, are obtained, as explained in the previous sec- 

tion, as minimal submodels from the faulty system model. A 
local state-parameter estimator is constructed for each iden- 
tification submodel and produces a local estimate 

p(x l f k , k \y l 0:k ) by using an appropriate algorithm. In this 
paper, we use an unscented Kalman filter (UKF) (Julier & 
Uhlmann, 2004) with a variance control algorithm (Daigle, 
Saha, & Goebel, 2012). 

4.2. Distributed Prognosis 

The local state -parameter estimates for each local distributed 
prediction module are constructed from the local estimates of 
the distributed fault identification submodels. Each predic- 
tion submodel is made up of a set of states Xi and parame- 
ters 0,, and constructs a local distribution p(x l j k ,9 l f k \y l 0:k ), 
by assuming that the local state-parameter estimates are suf- 
ficiently represented by a mean [i l and covariance E\ For 
each prediction submodel Mi(f ), we combine the estimates 
of the local identification submodels that estimate states and 
parameters in X r U0i into [i l and E\ If two submodels esti- 
mate the same state variable or parameter, then many different 
techniques can be applied depending on the desired perfor- 
mance of the prediction submodels, e.g., taking the estimate 
with the smallest variance, or taking an average. 

Several approaches can be used to perform prediction for each 


Algorithm 1 EOL Prediction 

Inputs: 

Outputs: 

for j = 1 to N do 

k <— kp 
*0) , i(j) 

x k x kp 

oT - oZ' 

while 7 % OI/ (xZ\ *l U) , ui) = 0 do 
Predict uj. 

Xfc+1 ~P(Xfc+lK W . 0 * O) ’Ufc) 

k <— k -f 1 
Y »0) , y *(j) 

x k x k + i 

n l ti) giU) 

U k ^ ^Jb+l 

end while 

EOL\ u) <- k 

Kp 

end for 


prediction submodel. In this work, given the mean and co- 
variance information, we represent the distribution with a set 
of sigma points derived using the unscented transform. Then, 
each sigma point is simulated forward to EOL, and we re- 
cover the statistics of the EOL distribution given by the sigma 
points (Daigle & Goebel, 2010). 

Algorithm 1 (Daigle, Saha, & Goebel, 2012), shows the 
pseudocode for the prediction procedure. The algorithm 
is executed for each submodel z, deriving local EOL pre- 
dictions using its local threshold function based on the lo- 
cal EOL constraints. For a given submodel, each sample 
j is propagated forward until T^ OLf (x t ^ k ,0^ k ) evaluates 
to 1. The algorithm hypothesizes future inputs uj.. Then, 
the global EOL/RUL is determined by the minimum of the 
local EOL/RUL distributions for each prediction submodel, 
i.e., p(EOL) kp |y j*, ) and p(RU L) kp |yj.* p ). To compute 
this, we sample from each local EOL distribution and take 
the minimum of the local samples. This is repeated many 
times and the statistics of the global EOL distribution are 
computed (Daigle, Bregon, & Roychoudhury, 2012). 

5. Case Study 

In this section, we apply our distributed diagnosis and prog- 
nosis approach to a four-wheeled rover testbed developed at 
NASA Ames Research Center. We develop a model of the 
rover, and demonstrate the approach using simulated scenar- 
ios. 

5.1. Nominal System Modeling 

The rover model was originally presented in (Balaban et al., 
2011). In this section we summarize the main features and 
include some extensions to the model. 

The rover consists of a symmetric rigid frame with four 
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independently-driven wheels. The wheel speeds are governed 
by 


VFL = "7 ( TmFL — Tjfl ~ r glFL + t 9 vFl) , 

JFL 

UFR ~ (jmFR -* TfFR ~ r glFR ~ T grFR.) , 

Jfr 

d>BL ~ — ( T mBL — TfBL “ r glBL + T grBL,) » 

Jbl 

UBR = {p’mBR - TfBR — TglFR ~ T grBFi) • 

Jbr 


(Cl) 

(C 2 ) 

(c 3 ) 

(C 4 > 


The F y By Ly and R subscripts stand for front y lefty back , 
and righty respectively. Here, for wheel w , J w denotes the 
wheel inertia; r mw is the motor torque; t/ w = [x j w uj w is 
the wheel friction torque, where fj,f w is a friction coefficient; 
T giw = r w {ji g i(v w — v) is the torque due to slippage, where 
r w is the wheel radius, fx g i is a friction coefficient, v w is the 
translational wheel velocity, and v is the translation velocity 
of the rover body; and r grw = r w fx grw u> cos 7 is the torque 
due to the rotational movement of the rover body, where n grw 
is a friction coefficient, u is the rotational velocity of the rover 
body, and 7 = arctan//6 with l being the rover’s length and 
b being its width. 

The translational velocity v of the rover is described by 

V = — ( FglFL + FglFR + FglBL + F q IBr) » (C 5 ) 


where m is the rover mass, and for wheel w y F g i w — fi gi (v w — 
v) is the force due to slippage. The rotational velocity u is 
described by 


w = — (d cos 7 Fg tF ft +dcosyF g iBR - dcosyF g i FL 

- d cos 7 FglBL - dF grFL - dF grFR - dF grB L 
~ dFgrBR )• (C6) 


Here, J is the rotational inertia of the rover and d is the dis- 
tance from the center of the rover to each wheel. 


The integral error terms are governed by 


e-iFL — ufl — ^FL» 

(Cll) 

eiFR = U FR — UFR, 

(C 12 > 

ZiBL = UBL ~ UBL, 

(C 13 ) 

CiBR ~ U BR - VBR- 

(C14) 


The batteries, which are connected in series, are described 
by a simple electrical circuit equivalent model that includes a 
large capacitance C in parallel with a resistance Rpy together 
in series with another resistance R s . 5 The battery charge vari- 
ables qi are governed by 


<h = —Vi/Rpi — ( i F L + IFR + 'BR + *Bl), (cis) 

<72 = - Vi/ Rj>2 - ( IFL + IFR + ^BR + *Bl)> (Cl6) 

<73 = — V 3 / Rp3 — {}FL + *FR + IBR + tBL), (C 17 ) 

<74 = -V^/RpA - ( iFL T i>FR T ^BR +»Bl)- (cis) 


The available sensors measure the voltages of the batteries, 


Vjf = qi/Cbi — Rsl * ( iFL + *FB + iBR + *Bl), (ci9> 

V 2 = Q2/Cb2 “ Rs2 * (^FL + *FB + tBR + *BI,)> (C 20 ) 

V 3 = <73/Cb3 - R s 3 * i^FL + iFR + *BH + *BI,)> (C2l) 

V 4 = q^/CbA - Rs 4 * (iFL + iFR + iBR + *Bl)i (C22) 

the motor currents, 

i*FL = iFL » (C23) 

i*FR ~ *FRi (C24> 

i BL— i BLi (C25> 

i*BR = iBR , (C26> 

and the wheel speeds, 

UJpL = UFL, (C27) 

U FR — U FR J (C28) 

U BL =UJ BL 1 ( c 29> 

uj BR = uj BR • (C30> 


The wheels are driven by DC motors with PI control that sets 
the voltages V applied to the motors. The motor currents i 

are governed by 


iFL = ~(VfL — ipL^FL ~ ku /<^FL)> 

( C 7 ) 

iFR = — (V>R - IFrRfR - fcu^Ffi)* 

(C8) 

iBL = ~[( V BL - iBLRBL ~ ku^BL)i 

(c 9 ) 

iBR = Jt(VbR ~ iBRRBR ~ k u u) B R), 

(Cio) 


where L is the motor inductance, 7 ? is the motor resistance, 
and is an energy transformation term. The motor torque 
is r Tnw = k T i W y where k r is an energy transformation gain. 
The voltages applied to the motors are determined by the con- 
trollers, where for wheel Wy V w = P * ( u ^ — lj w ) + / * e^, 
where P is a proportional gain, u w is the commanded wheel 
speed, I is an integral gain, and e{ W is the integral error term. 


Here, the * superscript indicates a measured value. 

5.2. Faulty System Modeling 

In this work, we consider different faults in the motors and the 
batteries. First, we consider friction-based damage progres- 
sion in the motors, resulting in an increase in motor friction 
over time. For wheel w y the fault progression function is de- 
fined as: 


A/FL = Vf FL VfFL V FL » 

(C 3 l) 

A/FB = VjpR VfFR V FR > 

(C32) 

ilfBL = VfBL V-fBL Wfli, 

(C33) 

2 

M/B/t = v jBR VfBR U B R, 

(C34) 


5 We use a simple model here only for demonstration purposes. More detailed 
battery models for prognosis can be found in the literature, e.g., (Saha & 
Goebel, 2009). 
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Submodel 

-~Xl 


T Ti 

Yi 

Ci 

Mi — 

9i 

Cbi 

^FLy^FR^BL^BR 

Vi 

Cl5,Ci9,C23,C24,C25,C2C 

M 2 

Q2 

c b2 

Vfl-> \ fr> i b L i i b r 

vi 

Cl 6 ,C20 >C23 ,C24 ,C25 ,C 2 6 

M 3 

93 

Cb3 

\fl > 2 fh> \bl > \br 

vi 

Cl 7 *C2 1 »C 2 3 ,C24 ,C 2 5 ,C26 

m 4 

94 

Cb 4 

^FL^FR^BL^BR • 

v: 

Cl 8 ,C 2 2 ,C23 ,C24 ,C 2 5 ,C26 

Ms 

iFL^iFL 

Rfl 

ufl,v fl 

*FL 

C7,Cn,C23,C27 

m 6 

i>FR, e-iFR 

Rfr 

UFRi^FR 

ip 

C8,Ci2,C24,C28 

Mi 

iBL,eiBL 

Rbl 

ubl^bl 

*SL 

C9,Ci3,C25,C29 

Ms 

lBR,ei B R 

Rbr 

UBRj^BR 

i*BR 

Cl0,Ci4,C26,C30 

m 9 

UFLyVyU, HfFL 

VfFL 

^FL^FR^^BL^BR 

u fl 

Cl,C3i,C5,C6,C23,C28.C29,C30 

M io 

UFR,V>U)yllf F R 

VfFR 

^FR^FL^BL^BR 

^fr 

C2,C32»C5,C6,C24,C27,C29,C30 

Mu 

UBL,V y U), fJifBL 

VfBL 

i B L > ^ p L i M f R i Wqr 

vbl 

C3,C33,C5,C6,C25,C27,C28,C30 

M 12 

VBR, VyUyflfBR 

VfBR 

^br^fl^fr^bl 

u br 

C4,C34,C5,C 6 ,C26>C27,C28,C 2 9 


Table 1 . Set of minimal submodels for the rover testbed computed from the nominal system model. 

Submodel Xj P, U x Y \ C\ 

Ms, 9 ^fTP^^T J^FLT lFLT^iFL ^/FL, Rfl ^FL , &FR » ^BL » U BR V FL > * FL U^UTJc) 

A^e.10 UFR,V,(jJ,HfFR,iFR,eiFR V/FR,Rfr UFR^FLi^BL^BR Wp R} i* FR CeUClO 

A't 7,11 FfBLiRBL U BL y<^FLi ^FRi^BR ^BL^BL C 7 U C\\ 

M&,12 ^BR,V,U),flfBRyiBR^iBR V/BR,RbR ^BR, V* fl , V fr , U>* bl U>bRi*BR Cq U C 12 

M 1,2,3, 4 91*92,93,94 Cbl,C b2 ,C63,C b 4 ^F L'^F R^ L *B L^ l *B R V\ ; ^ 2 * , ^3 iW C'l U C 2 U C 3 U C 4 


Table 2. Residual generation and fault isolation submodels. 


where /jbf w is the fault parameter, and Vf w is the fault pro- 
gression parameter. 

We also consider abrupt resistance increases in the motors, 
represented as an abrupt change in parameter R w for wheel 
w y with A R w as the fault progression parameter. 

For the batteries, we consider abrupt capacitance decreases, 
represented as an abrupt change in parameter Cbi for capacity 
i. A Cbi is the fault progression parameter. 

We are interested in predicting when any of the rover batteries 
are at their charge threshold, beyond which the batteries will 
be damaged. These faults can cause the charge thresholds to 
be reached earlier since they will affect current draw. The 
constraints are given as 


9l > 9 _ , 

(C 35 ) 

<12 > q ", 

(C36> 

93 > 9 _ , 

(c 37 ) 

94 >9", 

(C38> 

where the charge threshold is given by q~ 

= 2 x 10 4 C. The 

rover cannot be operated when any of the constraints C 35 -C 38 


are violated. 

5.3. Results 

To demonstrate the validity of the approach, we describe two 
different faulty scenarios of the rover. In the first, friction 
damage is progressing on one motor, and in the second, a 
capacitance decrease occurs in one battery. In all cases, the 
rover travels between various waypoints, moving at an aver- 
age speed of 0.5 m/s. Table 1 shows the minimal submodels 
for the rover derived by using measured values as local in- 
puts. Table 2 shows the submodels for residual generation 


Submodel 

~Xi 


~U~i 

t r 

~cc~ 

M 17 (^ 61 ) 

q\,Cb\ 

A<?61 

iFL,iFR,iBL, IBR 

0 

Cl5,Ci9,C35 

M 18 ( 0 , 2 ) 

92 , C b 2 

AC b2 

iFL,lFR , ii 3L, IBR 

0 

Cl6,C20,C36 

M 19 (^ 63 ) 

93, Cb3 

ACb3 

iFL,i F R,iBL, IBR 

0 

Cl7,C 2 l,C37 

•M2o(C64) 

94 , Cm 

A Cm 

iFL>lFR,iBL,iBR 

0 

Cl8,C22,C38 


Table 4. Prediction submodels for capacitance faults. 


and fault isolation. These submodels have been designed to 
obtain globally diagnosable subsystems by using the design 
algorithm in (Bregon et al., 201 1). In this work, we have con- 
sidered five subsystems, one for each wheel components and 
another one for the batteries. For example, the subsystem for 
the front left wheel components is not globally diagnosable if 
we only consider submodel Mq (which includes the front left 
wheel friction wear parameter, v/fl)- The design algorithm 
in (Bregon et al., 2011) determines that we need to merge 
submodels Mb and to make the front left wheel subsys- 
tem globally diagnosable. The process is similar for the rest 
of the subsystems. 

Minimal submodels for identification can be computed from 
the minimal submodels in Table 1 by defining the fault pro- 
gression function (if necessary), and by making the fault pa- 
rameter to become a state and the corresponding fault pro- 
gression parameter to become the parameter. Regarding pre- 
diction, the correct prediction submodels to use depend on the 
scenario, as will be shown later. 

5.3.1. Friction Damage Progression 

Wc first consider a scenario in which, for the front-left mo- 
tor, the friction begins to increase. The friction damage pro- 
gression begins at t = 50 s with friction wear parameter 
v fFL = 1 x 10 “ 3 s. A fault is detected by the local di- 
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Submodel 



TTi 

~y- 

~ci 

Ml3{^fFL) 

-Mis(m/fl) 

Mig^/fl) 

Ql , IFL , dFL, VFL > P'fFL 
<72 , IFL , ZiFL , C JFL , hfFL 
93, lFL,eiFL,U>FL,lJ'fFL 
q4,iFL,ZiFLiUFL, V'fFL 

Cb\,VfFL, Rfl 
C b2, v/fl,Rfl 
C b3,VfFL, Rfl 
C b4,VfFL, Rfl 

UFL,V , LJ, IFR.IBL , IBR 
UFL , V, U, IFR ,IBL,IBR 
UFL,V,(jJ, IFR, lBL,iBR 
UfL,V,UJ, IFR^BL , i BR 

0 

0 

0 

0 

Cl,C31,C7,Cn,Ci5,C35 

Cl ,C 31 ,C7,Cn ,Ci6,C36 
Cl ,C31 ,C7,Ci 1 ,Ci7,C37 
Cl,C 3 i,C7,Cn,Ci8,C38 


Table 3. Prediction submodels using commanded wheel speeds and rover velocities as local inputs. 




Figure 3. Current i* FL increase through time. 


agnoser computed from submodel M^g at 119.25 s, via an 
increase in the motor current i F i. The initial candidate list is 
immediately reduced to one candidate, based on the 

signatures and orderings (other faults in the front left wheel, 
like Rfl , produce different fault signatures). Thus the true 
fault is isolated. 

Fault identification is initiated once the candidate is isolated. 
For the friction damage progression fault, the wear rate is/fl 
estimate averages to v/fl = lx 10 -3 s with ver y small 
output error. Figure 2 shows the wear parameter estimate for 
friction damage. 

As a result of the continuously increasing friction, the current 
drawn by the motor increases as well in order for the motor 
controller to maintain the same desired wheel speed (Figure 
3 shows this increase in the current through time). Hence, the 
total current drawn from the batteries is increased, and EOL 
occurs around half an hour. Because ipi is constantly chang- 
ing, and in a way that is dependent on the motor state, it is 


incorrect to use it as a local input for prediction and to de- 
compose the prediction problem into independent local pre- 
diction problems for the batteries and motors, i.e., it is not 
known a priori. Therefore, we compute submodels using as 
local inputs average values for the remaining motor currents, 
average commanded wheel speeds, and average rover trans- 
lational velocity v and rotational velocity uj. The prediction 
submodels for this case are shown in Table 3. EOL for this 
fault is computed by merging the local EOL from those sub- 
models in the table. Note that the prediction submodels used 
in this case do not correspond directly to those used for es- 
timation. So, when constructing the estimate for .M13, for 
example, it takes the estimates from M\ and Mg. 

The prediction results are shown in Figure 4. The increased 
friction causes the batteries to discharge faster, and EOL oc- 
curs around 1650 s. Here, we used the relative accuracy (RA) 
as a measure of prediction accuracy, and the relative standard 
deviation (RSD) as a measure of spread. Each prediction met- 
ric is averaged over multiple prediction points (one every 100 
s of usage) (see (Saxena et al., 2010; Daigle, Saha, & Goebel, 
2012) for the mathematical definitions of these metrics). For 
this experiment, RA averages to 91.63% and RSD averages 
to 16.26%. 

For the sake of comparison, we also ran this experiment us- 
ing the centralized approach. Figure 5 shows the prediction 
results obtained. Looking at the prediction metrics, we see 
that the centralized approach behaved very similar to the dis- 
tributed approach but a little bit worse, with RA averaging 
90.90% and RSD averaging 17.72%. However, this is just 
a particular example, but, in general, both approaches obtain 
equivalent results. 

5.3.2. Capacitance Decrease 

As a second scenario, we consider a capacitance decrease 
fault in battery 3 of the rover, C 03 . The fault begins at t = 50 
s with an abrupt decrease from 2000 to 1800 in the capacity 
of the battery. The fault is detected immediately by the local 
diagnoser computed from submodel M 1,2, 3, 4 at 50.0 s, via an 
increase in the voltage V3. The fault candidate is immediately 
isolated, {C63}, based on the signatures and orderings, thus 
starting the fault identification. For the capacitance fault, the 
estimated value of the capacitance averaged Cb 3 = 1798.6 
C with very small output error. As a result of the decrease 
in capacitance, the battery discharges at a faster rate, and so 
reaches end of discharge more quickly. The prediction sub- 
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Figure 4. Predicted RUL of the rover for the distributed ap- 
proach. The mean is indicated with a dot and confidence in- 
tervals for 5% and 95% by lines. The gray cone depicts an 
accuracy requirement of 15%. 
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Figure 5. Predicted RUL of the rover for the centralized ap- 
proach. The mean is indicated with a dot and confidence in- 
tervals for 5% and 95% by lines. The gray cone depicts an 
accuracy requirement of 15%. 


models for faults in the capacity of the batteries are shown 
in Table 4. For this scenario, with a fault in C& 3 , we used 
submodel obtaining RA average to 98.25% and RSD 
average to 10.12%. 

6. Conclusions 

This paper presented a distributed integrated model-based 
diagnosis and prognosis framework. Our approach starts 
off with a common modeling paradigm to model both the 
nominal behavior and fault progression, and then proposes 
a framework where the global system model is decomposed 
into smaller independent submodels. These submodels are 
then used to distribute the different diagnosis and prognosis 
tasks. Model decompostion is carried out based on the re- 
quirements and constraints of each task. We demonstrated 
our approach on a four-wheeled rover testbed, where we di- 
agnosed faults and prognosed the EOL/RUL accurately. We 
compared results obtained by using our distributed approach 
against those obtained using a centralized approach, showing 
that both approaches obtain the same results. 

Most approaches in the literature focus in either the diagnosis 
or the prognosis task. Some works have proposed the inte- 
gration of both tasks within a common framework (Patrick 
et al., 2007; Orchard & Vachtsevanos, 2009; Roychoudhury 
& Daigle, 2011), however, unlike our approach, these ap- 
proaches perform the diagnosis and prognosis tasks in a cen- 
tralized way, thus suffering from scalability issues due to the 
large number of states and parameters in real-world systems. 
To the best of our knowledge, there is no approach in the lit- 
erature which combines, in a distributed way, the diagnosis 
and prognosis tasks. Our approach is limited by the number 
and location of the sensors in the system. Since our decom- 
position algorithm is guided by the set of available sensors, 
the distribution capabilities of the approach is determined by 
them. 

In future, we will apply this approach to larger systems, to 
study the scalability of our diagnosis and prognosis scheme, 
and will perform a more detailed comparison against the re- 
sults obtained by using a centralized approach. We will also 
extend the capability of this approach to hybrid systems, as 
well as diagnosis and prognosis of multiple faults. 
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